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1. Introduction 

In this work wc deal with the standing wave solutions 

of the time-dependent Schrodinger equation with electromagnetic field 
= (f V - A{x))\ + W{x)^ - Iv^r V, 
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where the Schrodinger operator is defined as 

(|V-yl)^ := ~H'^A~f (A I V) + l^ll^ - f divA 

Here > is the Planck constant, p G (1,5), and the functions VF: K"^ ^ M and 
A: — > are, respectively, a scalar potential of the electric field E = — VW and a 
vector potential for the external electromagnetic field B = curl A. Now, the function 
w: R'^ C which appears in (p{x,t) satisfies, more generally, a time-independent 
equation of the form 

(|V - A{x))^u + V{x)u = K{x)f{\u\^)u, (1.1) 

where V{x) = W{x) + Vq, K : M.^ ^ R is an additional potential function, and 
/ : R+ — > R is a suitable nonlinearity. Quite recently, under reasonable assumptions 
on A, V and K, the study of the existence of ground (bound) state solutions un 
to (1.1) and the related investigation of the semi-classical limit (the transition from 
Quantum to Classical Mechanics as 7i — > 0), has been tackled in various contri- 
butions (see e.g. [2,4,6,7,14,19] for the case ^ 7^ and [3,11,12,13,15,23,26] for the 
case ^ = 0). More precisely, it turns out that, if zo € R"^ is a non-degenerate critical 
point of the so called ground-energy function : R'^ ^ R (see Definition 2.2), then 
for every H sufficiently small (1.1) admits a least energy solution uh concentrating 
near zq. In the opposite direction, we are interested in discussing some necessary 
conditions for the concentration of a sequence of bound-state solutions to (1.1) in 
the neighborhood of a given point zq- In absence of the electromagnetic field, this 
problem has been studied in various papers (see e.g. [1,29,30]), mainly in the case 
where f{u) is a power of exponent p (see also [17,24]). It turns out that, at least 
in this particular situation, for the concentration to occur, zq has to be a critical 
point for the ground-energy map (see [30, Lemma 2.5]) 

E^(z) = Zil-lM, for every e R^ (1.2) 
Kp-^ (z) 

On the other hand, to our knowledge, for a more general nonlinearity f{u), the 
function S,.(-) is locally Lipschitz continuous, and its further smoothness properties 
seem to depend on the uniqueness results for the limiting equation 

-Au + V{z)u = K{z)f{\u\'^)u, (1.3) 

where z G R^ acts as a parameter. To overcome this problem, recently, the authors 
have provided in [28] new necessary conditions involving generalized derivatives of 
Sr such as the Clarke subdifferential or even weaker conditions, not requiring any 
regularity of S^. (see Definition 2.4). 

Our purpose in this paper is to understand what happens under the presence 
of an external electromagnetic vector potential A, and to see whether A may in- 
fluence or not the location of spikes for the solutions of (1.1). Actually, in general, 
this fact seems to depend on the notion of concentration that one adopts. We con- 
sider at least two ways of saying that a sequence (u^) of bound-state solutions 
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to (1.1) is peaking around a given point zq. The first one, the most intuitive, is 
a pointwise concentration and it is precisely the one used in two papers by Wang 
and Zheng [29,30]. The second is a sort of energetic concentration in terms of the 
functional associated with (1.1), 



As we prove in the main result, Theorem 3.1, the vector potential A might affect 
the location of pointwise concentration points, whereas it does not influence the 
energetic concentration points. In the particular but fairly significant case where / 
is a power nonlinearity, the above notions of concentration coincide (see Proposi- 
tion 2.1), and it turns out that the peaks locate at the classical critical points of 
the smooth function (1.2) independent of A, thus rigorously confirming what con- 
jectured in [7] . In some sense, from an heuristic point of view, A tends to lurk into 
the complex phase factor of the solutions. We point out that, in the course of the 
proof of Theorem 3.1, we will derive an ad-hoc Pucci-Serrin type identity for the 
complex-valued solutions to (1.1) (cf. formula (3.6)). Just for the sake of simplicity, 
we restrict the attention to the physically relevant case of space-dimension n = 3. 



(1) 5Rw (resp. 3w) stands for the real (rcsp. the imaginary) part of w S C. 

(2) i is the imaginary unit, namely i^ = —1. For w € C, we set w = — iSw. 

(3) The gradient of a function / : — > M will be denoted by V/. The jacobian 
matrix of a function g : — > will be indicated by g'. The directional 
derivatives of / and g along a vector w will be indicated by ^ and ^ . 

(4) {x I y) denotes the standard scalar product in of x and y. 

2. Problem Setting and Auxiliary Results 

In this section, we collect a few preliminary definitions and results that we need 
in order to state and prove our main achievement. Theorem 3.1. For the sake of 
simplicity, we rename the constant H into £ > 0. We assume that the functions 

yl:R^^R^, F:M^^IR, JfiM^^M 

are all of class C^, K is positive and there exist Vo > and Kq > with 



Jf^{u) = i / \D^uf + V{x)\ufdx - [ K{x)F{\u\'^)dx, 




where = i 



V — A[x). Precisely, wc require that 



lim h ^Jfi{un) = Sr(-2o)- 



Notations 



inf V(x) = Vo 



and 



sup K{x) = Kq. 



(2.1) 
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Moreover, the function / : IR+ ^ R is of class C^, increasing, /(O) = and 
f(s) 

lim = and < i?F(s) < f{s)s for some p e (1, 5) and § > 2, 

where F{s) = ^ f{t)dt for s € M"*". In order to formuhitc the problem in a 
suitable variational setting, for every e > 0, we introduce the (real) Hilbert space 
y defined as the closure of C^(IR^, C) with respect to scalar product 

{u,v)n%y-=^ D^uWv + V{x)uvdx, D^u = ^V-A{x). 

As remarked in [14], y has in general no relationships with H^{9?, C). However, 
the following diamagnetic inequality is well known (see e.g. [21]) 

£:|V|u|(a;)| < \D^u{x)\, for every u 

gHav and a.e. xeR^, (2.2) 

so that juj e -ff^(M^,M) for any u G y. Finally we recall that the Schrodinger 

operator is gauge invariant: if we replace Ahy A = A + Vx for any x S C^(M^, R), 
and we let ii — ei^u, then curl A — curl A and 

(f V - i) w = ei^ (f V - A) u, 

so that IImII^^--^ ^= 

Under the above assumptions, we give the following 

Definition 2.1. We say that (u^) is a sequence of bound-state solutions to 

(f V - A{x))\ + V{x)u = K{x)f{\u\')u {S,) 
if Us belongs to 'H'% y for every s > 0, 

supe-^jjuejj^e^^ < oo (2.3) 

and Ue satisfies {Sg) on in weak sense. 

2.1. The ground- energy functions 
Fixed z &M.^, we consider the functional 

Iz{u) = \l |Vw|2 + V{z)\u\^dx - / K{z)F{\u\^)dx 

associated with the limiting equation (1.3). It is readily seen that Iz is over both 
the spaces //^(M^^M) and iJi(M3,C). 

Definition 2.2. We define the real and the complex ground-state functions 

: ^ M and Sc : ^ M 
by setting, for every ^ e M^, 

T,r{z) = min Iz{v) and ^c{z) = min Iz{v), 
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where Afz (resp. Afz) arc the real (resp. the complex) Nchari manifolds, 

A/; = {u e i^^(R^R) \ {0} : I'JuM = 0} , 
M, = {ueH^{R\C)\{Q}: I'z{u)[u] = 0} . 

Here li, (u) [v] stands for the directional derivative of Iz at u along v. 

We denote by Sr{z) the set of positive radial solutions up to translations to (1.3) 
at the energy level T,r{z). As the next lemma claims, the map enjoys some useful 
regularity properties (see [30]). 

Lemma 2.1. The following facts hold: 

(i) Sr is locally Lipschitz continuous; 

(ii) the directional derivatives from the left and the right ofY^r at every point z gM.^ 
along any w gR^ exist and it holds 

(^] W= sup {Vzlz{v)\w), 
\ aw J ves,.{z) 

Explicitly, we have 

\dwj veSr{z)^dw Jr3 2 dw J^s J 

{^] iz)= inf — (z) / ^(ix-— (z) / F(\vf)dx , 
for every z,w e . 

The next result will turn out to be pretty useful along the proof of our main 
theorem. We stress that it contains, as a particular case. Lemma 7 of [19]. 

Lemma 2.2. The following facts hold: 

(i) Y.ciz) = T^riz), for every z €M.^; 

(ii) if Uz : M.^ ^ C is a least energy solution of problem (1.3), then 

\V\Uz\{x)\ = |V;7^(a;)| and ^iUz{x)VUz{x)) = 0, 
for a.e. x G M^; 

(ill) there exist u) gR and a real least energy solution Uz of problem (1.3) with 

Uz{x) = e"^Uz{x), for a.e. x G (2.4) 
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Proof. Fix z e M'^. For the sake of convenience, we introduce the functionals 
T{u) = I \Vu\'^dx, 

P,{u) = ( \k{z)F{\u\'') - \V{z)\uf\dx. 

Observe that Iz{u) = \T{u) — Pz{u). Consider the following minimization problems 

ar{z) = min{r(?i) : u e i^^(R^M), Pz(u) = l}, 
ac{z) = min{r(7i) : ?i G i^l(R'^C), Pz{u) = l}. 

Note that, obviously, there holds (Jc{z) < ar{z). If we denote by the Schwarz 
symmetric rearrangement (see e.g. [3,21]) of the positive real valued function \u\ G 
H^{M.^,R), then, Cavalieri's principle yields 

/ F{\u^f)dx ^ / F{\uf)dx and / \Ui,fdx = / \u\'^dx, 

JR3 JR3 JR3 

which entails Pz{ui,) — Pz{u). Moreover, by the Polya-Szcgo inequality, we have 

T{u^)= [ |Vw*pda;< / |V|u||2dx < / \Vu\'^dx = T{u), 
Jk^ Jr'^ Jr3 

where the second inequality follows by (2.2) with A ~ and e = 1. Therefore, one 
can compute adz) by minimizing over the subclass of positive, radially symmetric 
and radially decreasing functions u & H^{M.^,R). As a consequence, ar{z) < <Jc{z). 
In conclusion, ar{z) = <Tc{z). Observe now that 

T,r{z) = min{7;,(u) : u e iJ^(M^,IR) \ {0} is a solution to (1.3)}, 
Y.c{z) = min{7^(u) : u e iJ^(M^,C) \ {0} is a solution to (1.3)}. 

The above equations hold since any nontrivial real (resp. complex) solution of (1.3) 

belongs to Mz (resp. A4) and, conversely, any solution of T,r{z) (resp. 'Sc{z)) pro- 
duces a nontrivial solution of (1.3). Moreover, it follows from an easy adaptation 
of [3, Theorem 3, p. 331] that T,r{z) = <Jr{z) as well as T,c{z) — adz). In conclusion, 

T,r{z) = ar{z) =^ adz) = 5]c(z), 

which proves (i). To prove (ii), let Uz : R^ ^ C be a least energy solution to 
problem (1.3). There holds |V|f/0|| < |Vf/0|. Assume by contradiction that 

C^{{x €M.^ : \V\Uz\ix)\ < \VUz{x)\}) > 0, 

where is the Lebesgue measure in K^. Then we get Pz{\Uz\) = PziUz) and 

ar{z)< [ \V\Uz\\^dx< [ \\/Uz\^dx = adz), 

which is a contradiction. The second assertion in (ii) follows by a direct computa- 
tion. Indeed, a.e. in M'', wo have 

|V|J7^|| = \yUz\ if and only if mzV{QUz) = QUzV{mz). 
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If this last condition holds, in turn, a.e. in we have 

which implies the desired assertion. Finally, the representation formula of (iii) is an 
immediate consequence of (ii), since one obtains Uz — e^'^lC/^j for some w e M. □ 

2.2. Generalized gradients 

Assume that / : — > M is a locally Lipschitz continuous function. For the reader 
convenience, we recall that the Clarke subdifferential (or generalized gradient) of / 
at a point z (cf. [8] ) is defined as 

dcf{z) = I?? e : f{z, w)>{r]\w), for every w e M^j, 

where f'^{z,w) is the Clarke derivative of / at along the direction w, defined as 

J (z; w) = limsup ^ —. 

Prom [8, Proposition 2.3.1] we learn that dcf{z) is nonempty, convex and 

dci~.f){z)^~dc.nz), for every ^ e (2.5) 
In light of (i) in Lemma 2.1, we arc allowed to give the following 

Definition 2.3. We denote by 6 C M"^ the set of critical points of the function 
in the sense of the Clarke subdifferential, namely 

S := {z G : € dcT^riz)}. 

Now, for z € M^, we consider the gauge invariant functional Jz : H^{R^, C) — > K 

J,(u) = i/ \{lV-A{z))u\' + V{z)\u\''dx-[ K{z)F{\u\'')dx, 

associated with the limiting equation 

(iV - A{z)f u + V{z)u = K{z)f{\u\'^)u. 

We denote by Gc{z) the set of the nontrivial solutions : M"^ ^ C, up to transla- 
tions, of the above limiting problem with bounded, but not necessarily least, energy. 
Moreover, we introduce the linear map : — > M, defined as 

3 

"^zix) := Aj(z)xj, for every a; S M"^. 

Apparently, for every ^; e K^, there holds z{x) = A{z). It is readily seen that 
for every v G Gc{z) wc can write v = e^'^'Uz, where Uz is a (possibly complex) 
solution to problem (1.3). 
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Definition 2.4. Let z e K.^. For every w E M.'^ we define (w) and Tf{w) by 

F;H:= sup {y,Mv)\w) and F+H := inf {V,J,{v)\w), 
veGciz) veGc{z) 

wliere Vz is tlie gradient with respect to z. Explicitly, for every w gM.^, 

r / OA f - \ 
T-{w)= sup {jr^z) \ / ^{iU,VU,)dx) 

r+H = - inf [(I^Wl/ ^(ilJz^Vz)dx) 

dw J^3 2 Sw; 7ii3 ' ' . 

Notice that 

dTf{0) = |?7 e R3 : F^(w) > (r/ | w;), for every w e M^j, 

where 9F^(0) is the sub differential of the convex function F^ at zero. We set 

S* := e : e dT-{0) D 9F+(0)} 
and we say that S* is the set of weak-concentration points for problem {Ss). 

2.3. Concentration of bound-state solutions 

We now introduce two (gauge invariant) notions of concentration for a sequence of 
bound-states solutions of (S^) around a given point. 

Definition 2.5. Let zq e R'^ and assume that (we^) C H^Ay is a sequence of 
bound-state solutions to problem (Se). We say that 

(i) Zq is a concentration point for (we^) if |ue^(2;o)| > P > and for every > 
there exist p > and /lo > 1 such that 

|'Ue^(x)| < rj, for every h> ho and |a; — zq] > ShP- 

The set of such points will be denoted by C M^; 

(ii) zo is an energy- concentration point if 

lim s^^Je^iu^^) ^ Y.r{zo). 

n — ^oo 

The set of such points will be denoted by # c M^. 

For instance, if K = 1, / is a power, zq is a minimum point of V and (we^) is a 
sequence of least-energy solutions to (Se), then zq S ^$ (cf. [19, Lemma 3]). 
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Next we see that in the case of power nonhnearities 

f(u) = for some p G (1,5) and A > 0, (2.6) 

the above notions (i) and (ii) coincide. 

Proposition 2.1. Let f be as in (2.6). Then S = . 

Proof. Let zq ^ S and consider Vh(x) = Ue,^{zQ + Shx). Then {\vh\) converges to 

some w > weakly in ^^(IR^, K) and strongly in L\^J^ ■, R) for 2 < g < 6 (see Step 
I in the proof of Theorem 3.1). By Kato's inequality [27, Theorem X.33], we get 

/ K(zQ-\-ehX)\vh^vdx> \ V\vh^v -\-V{zq+ £hx)\vh\vdx 
Jr^ Jr3 

which, as /i — > oo, yields, 

/ K{zoMP+'^dx> [ \Vv\^ + V{zo)\i\'^dx. 
Therefore, there exists i? € (0, 1] such that € Afz^- As a consequence, 

Sr(^o) < (i - ^) IVSI^ + V{zo)\vfdx 

^ (5-^)lrinf / \V\vh\\^+V{zo+ehx)\vh\^dx 

^(h- ^) lipiiif / I(t^ - A{zo+ehx)) Vhf + V{zo + Shx)\vhfdx 



< liminfe^^JeJweJ = ^^(^o), 

n — i-oo 



where we have used the diamagnetic inequality (2.2) with e = 1. Hence we get 
■!? = 1, which gives at once v G Afzg- Then, 

/ + V{zoM'dx < U - ^y\imMe^'j,,M 



(i - ^)"'s.(^o) < \^vf + V{z^)\v\^dx. 



This implies that \vh\ v strongly in iI^(IR^,IR). Repeating the arguments in 
the proof of [19, Lemma 5] we conclude that zq G ^ (the concentration occurs 
exponentially fast, see Step II of the proof of Theorem 3.1). This proves that S (Z^ . 
The converse inclusion follows by the uniqueness of solutions (up to translations) 
to problem (1.3). Indeed, if Zq G ^, the sequence e^^Je^(M£^) converges to Jzo{vo) 
being vo an element of the family 

{e^^^o(-)+-</,o(a;)}^eM, 

where is the unique solution to (1.3) up to translations (cf. [19, Lemma 7]). In 

particular, there holds Jza{vo) = Izoi't'o) = ^r{za), that is zq G S", concluding the 
proof. For similar considerations in the case A = 0, see e.g. Lemma 4.2 in [18]. □ 
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We are naturally lead to consider the following question (see also Remark 3.1). 
Question 2.1. When f{u) does not satisfy (2.6), is it still true that (o = 

3. The Main Result 

For every p e (1, 5), let us set 

6p := {z e : {5 - p)K{z)\7V{z) = W{z)VK{z)}. 
We now come to the main result of the paper. 

Theorem 3.1. Assume that there exist C > and 7 > such that, for \x\ large, 
\A'{x)\ < Ce^l^l, |VV(a;)| < Ce'>'l^l, \VK{x)\ < Ce^l^L (3.1) 
Let (ue^) C H^J'y be a sequence of bound-state solutions to {Se)- Then, 

<^ C 6* and C &. 

If in addition f satisfies (2.6), then we have 

= ,r C 6 = 6p. 

Proof. Let zq £ 'rf and set Vh{x) = Uef^{zo + Shx) for every h > 1 and x e . 
Then, the sequence {vh) satisfies the rescaled equation 

-Avh - J {A{zo + Shx) I Vvh) - ^div^(2;o + ehx)vh 

+ \A{zo + ehx)\^Vh + V{zo + ehx)vh = K{zo + ehx)f{\vhf)vh. (3.2) 

We shall divide the proof into five steps. 

Step I. Up to a subsequence, {vh) converges in some Holder space Cf^"{M.^) to the 
function Vo{x) = e^'^^o^^^Uzoix), where Uzq : ^ C is a solution to the equation 

-AUz, + V{zo)Uz, = K{zo)f{\Uz,\^)Uz,. (3.3) 

By the assumption on (ug^), the sequence (vh) is bounded in H\ y, and the dia- 
magnetic inequality (2.2) immediately implies that {\vh\) is bounded in i?^(R'^,IR). 
Therefore, up to a subsequence, it converges weakly in i?^(IR^,R) and locally 
strongly in any L'(K^,M) with g' < 6 towards a positive function u*. Moreover, 

for each compact subset A C R'^, by the continuity of A, (vh) is also bounded in 
H^{A, C). We may now use the subsolution estimate (see e.g. [16, Theorem 8.17]) to 
get that {vfi) is also bounded in LJ^p(R^) and hence in Cf^"{M.^), via Schauder' esti- 
mates. By combining this fact with the results of [20], up to a subsequence, Vh con- 
verges to vq in Cj'^'"(M^) and furthermore vq ^ 0, since |?Vi(0)| = |u£^(2o)| > £• > 0. 
By continuity, the limit vq satisfies the limiting equation 

-Avo - f {A{zo) I Vvo) + \Aizo)\^vo + V{zo)va = K{zo)fi\vo\'')vo. (3.4) 
If we define : a; € ^ e-''^-o(^^vo{x), then U^o satisfies (3.3). 
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Step II. There exist two positive constants i?* and C* such that 

|fh(a;)| < C*e-V^l^l, for every |x| > i?* and /i > 1, (3.5) 

where Vq is defined in (2.1). Since zq G "^^j we have Wft(x) as |x| — > oo, uniformly 
with respect to ft, > 1. Hence, for any > 0, we can find a radius i?^ > such 
that < t] whenever \x\ > Rr^ and h>l. Therefore, exploiting again Kato's 

inequality 

A|t;/i| > ^{vh\vh\~^{S ~ i^)^Vh) (in distributional sense), 
and taking into account that / is increasing, there holds 

A\vh\>V{zo + ehx)\vh\-K{zo + ehx)f{\vhf)\vh\ > [Vo - Kof{v'')]\vh\ 

in the sense of distributions on {|x| > Rr,}, where ii'o > is as in (2.1). Let 

To be a fundamental solution for —A + c,,, where c,, = Vo — Kof{rf'). We can 
choose To so that |w/i(x)| < [Vq — KQf(Tf )]TQ{x) holds for = i?^. Then, if 
w = \vh\ - [Vo - KofiTj^)]To, there holds 

Am; = A|z;;,| - [V^ - ifo/(?7')] AFq 

> [Vo - Kof{v^)]\vh\ - [Vo - Kof{v^)fro 
= [Vo - Kof{rf)]w 

in distributional sense over {|a;| > R.q}. Then, by the maximum principle, w{x) < 
for every \x\ > Rn- Since, as known, Fq decays exponentially at the rate y/c^, fixing 
1] = 7]^ so small that /(r;^) < Vq/2Ko, we can find constants i?, > and c > such 
that ro(x) < cexp{— A/Vb/2|a;|} for \x\ > R*, which yields the desired conclusion. 
Step III. For every ft > 1, the following identity holds 

(^(^o + e/ia;) | A{zo + ehx)J\vh\'^ - ?li(^jVvh \ -^{zo + ehx)vhj 
— {zo + ehxy-^-—{zo + ehx)F{\vh\^) 



I 



dx = 0. (3.6) 



Rigorously, we cannot directly apply the Pucci-Serrin variational identity [25], since 
the solutions to equation (3.2) are complex-valued. For we are not aware of any 
explicit reference to cite for the identity we need, we will derive (3.6) directly (see 
also [9] ) . Throughout the rest of this step only, we use the less cumbersome notation 
a; • y in place of {x \ y) to indicate the standard scalar product in 9?. 
First of all, let us observe that, for every ft > 1, 

|Vt;,,| < \D^Vh\ + \A{zo + ehx)\\vh\. 

Hence, taking into account Step H and the bounds (3.1) and (2.3), we get 

||Vwhj|L2(R3) < \\D^Vh\\L^(w) + \\A{zo + ehx)vh\\L'^(m) 

< WD^vhWl^^^) + \A{zo)\\\vh\\LHv.^) + c\\e-''^\^^\x\vh\\LHvJ^) < c, (3.7) 
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for all /i > 1 and some c > 0. Let S > and consider the cut-off function ipg = ■4>{6x), 
where tjj e Cc(R^) is such that V(a;) = 1 for \x\ < 1 and ijj{x) = for |a;| > 2. If 

gfe denotes the fc-th vector of the canonical base in R^, we test equation (3.2) with 
the function ipsSk ■ and we take the real part. Firstly, we have 

^ [ Vu/, • V[V'5efc • Vvh]dx = / Vvh ■ ^ipsek ■ Vvhdx - f Vip^ ■ ej-^^^^^dx. 

Jrs Jrs ^ 

As a consequence, by virtue of (3.7), the Dominated Convergence Theorem yields 
lim 3? / Vvh ■ V[V'5efe • V^dx = 0. 

Now, we have 

3? / K{z,i + ehx)]{\vh'^)vhi\)be.k ■ Vt^rfa; = 5i / K{zq + ShxJ^sek ■ VF{\vhf)dx 

= -£h[ ^{zo + ehx)^sF{\vh\'')dx- [ K{zo + ehx)^F{\vh\^)dx. 
V dxk V 'dxk 

Hence, in light of (3.1), (3.5) and (3.7), by the Dominated Convergence Theorem 

we have 



limK / K{zQ+ehx)f{\vh\^)vhil^5ek-Vvhdx = -£h [ ^^{zo+ £hx)F{\vh\^)dx. 
In a similar fashion, there hold 

lim / V{zq + £hx)vh'4)&ek ■ Vvhdx = -Sh / ^ — (-zo + £hx)—^dx, 

lim 3?/ \A{zQ+ehx)\'^Vhil^5ek-Vvhdx = ~£h I A{zQ+£hx)-- — {za+£hx)\vh\^dx. 
Finally, we have 

J{5) = 3? / f A(^o + Shx) ■ ^Vhi)sek ■ W^dx = Ji{d) + J2{5) + JaiS), 

JR3 



where we have set 



Jl{S) = -£h^2^ I j-g^{zo+£hx)tps-g^Vhdx, 



m=l 

3 



J2(5) = -3? V / f A„(0o +e/ta;)^-|^?;hdx, 
/m3 ^Xto 



J3(5) = -3? > ' / jA^{zo + £hx)'ips^ — ^^Vhdx. 



3 

/r3 ' 5xfe9a; 
After a few computations, one shows that J2(i5) ^ as ^ ^ and 



MS) = -3? / 2£!LdivA(zo + £hx)vh^5ek ■ Vvhdx - J(<5) + Q{5), 
Jm3 
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with Q{S) ^ as 5 — > 0. Furthermore, again by (3.1), (3.5) and (3.7) 

f OA 
hm Ji{5) = -£^SR / jVvh ■ ^ — (^^o + £hx)vhdx. 
5^0 Jrs ' dXk 

Therefore, we obtain 

UmJ(^) = -K/ ^diyA{zo+ehx)vhek-Vvhdx-Sh^ T^^ft""^ — {zo+£hx)vhdx. 

J^3 Jr3 OXk 

Adding the above identities immediately yields (3.6). 

Step IV. We apply the Dominated Convergence Theorem to take the limit as 
/i — > oo into identity (3.6). The only troublesome term is 

^(^j'^Vh I ^^{zo + ehx)vhy 

since wc apparently have no control on the decay of Vu/j. Taking into account (3.7) 
and recalling that \/vh{x) Vvo{x) for all a; G R"^, up to a subsequence, we have 

Wvh^'Vvo, weakly in i^(M^). (3.8) 

On the other hand, by virtue of Step II, there exist i?* > and c > such that 

< ce-(^ ^-7^01^1^ for every |a;| > i?*. 



dA 

g^izo + ehx)vh 



Consequently, since Vh{x) — > vo{x) for all a; e and A e C^(]R^), there holds 
dA dA 

-^{zo + ehx)vh ^ -^{zq)vo, strongly in L2(IR^). (3.9) 

Thus, by combining (3.8) and (3.9), for each fc, we immediately get 

f /I dA \ [ I dA \ 

lim / '^i-\/vh\-^ — {zo + ehx)vh)dx= ^(j'^'Vo \ ^ — {zo)vo)dx. 
Jr3 \i dxk I Jw ^ ' dxk / 

Since similar considerations apply to the other terms that appear in (3.6), we can 
therefore pass to the limit as /i ^ oo, to find, for each k, 

I dA \ I dA \ 

(^{zo) I A{zo))\vo\^ - ^\Vvo I ^{zo)vo)dx 

]j3 \ OXk I \ OXk ' 

dV , . f |vopj dK 



/ ^-dx-^{z^)! F{\vo\'')dx = Q. (3.10) 

JR3 ^ OXk Jr3 



Now, as proved in Step I, vq can be represented as vo{x) = e*^^o C/jq (a;) where 
: — > C solves (3.3). Taking into account that 

iVvoix) = e'^'oi-)A{zo)U,,{x) - ie'^'oi-)yu,,{x), 
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for every x G K'' we obtain 
dA 



(^(^o) I A{zo))\vo{x)\' - ?fi{lVvo{x) I §^^{zo)vo{x) 
= {^^izo)\A{zo))\U,dx)\' 



dA 
\dxk 



= /^{zo)\A{zo))\U,,{x)f 



((^M I Aizo))\U.,{x)\' - ^{^^{zo) | VC/.„(a;))i7.„(x) 
dA 



./ dA 



= {^^{zo)\WU,,{x)VU,,ix)) 

Hence, equation (3.10) can be rephrased as 
/ dA 



\dxk 



(^o)l 



^R{iU,,^U,,)dx 
dV 



2 dxh 



(zo) / i^(|C/.j')da; = 0, 



/ dA 

I —{zo) I / 3?(iC/.„VC/,„)dx 



for every A; = 1, 2, 3, namely, 

idA 

\dw 
for every w e K^. 

Step V. In this final step, we prove the desired inclusions stated by the theorem. 

As a consequence of identity (3.11), in light of the definition of r^{zo;'w), we im- 
mediately deduce that zq E &* , thus proving that C 6*. Let us now assume that 
zq € (a. Then Jzo{vo) = ^c{zo) = ^r{zo), and by virtue of (iii) of Lemma 2.2, we 
have Uzo{x) = e'^'^Uzg{x) for some w e K, where u^g is a real least energy solution 
to (1.3). Moreover, by (ii) of Lemma 2.2, we have 

5i {iUza{x)'^Uza{x)) = 0, for a.e. x e M^. 

Then, in light of Lemma 2.1 and (3.11), we obtain 



as 



sup 



. \u\^ dK 



= sup 
dV 



dV_ 
dw 

dV 
dw 



F{\u\^)dx 

(zo) [ \^dx+^{zo) I F{\U\^)d: 

( \E^dx+^{zo) [ Fi\Uz,\')dx = 0, 
Jr3 I aw J^3 
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for every w e R'^. In a similar fashion, there holds 

for every w e M^. In particular, by the definition of {—'Er)^{zo; w), we get 

(-i],)°(^o;H> (^^)Vo)>0, 

for w e M^. Hence e dc{—'Sr){zo), which, in light of Proposition 2.5, yields 
Zo € 6. Finally, if f{u) satisfies (2.6), problem (1.3) admits a unique real solution 
(f)o up to translations (see [5]). Taking into account Lemma 2.2, there exists w G R 
such that vo = e'^-o(^)+''^(/)o(a;). Then, if e = ^ (see Proposition 2.1), we 
have <S'r(^o) = {^o}, admits all the directional derivatives and, by the above 
inequalities, 

for w G R^. Since up to a multiplicative constant S^. writes down explicitly as (1.2), 
the last assertion readily follows by a direct computation. □ 

In light of identity (3.11), we also have the following 

Corollary 3.1. Under the assumptions of Theorem 3.1, for every zo S there 
exist constants Ai,A2, A3 G M (possibly zero) and 71,72 € M\ {0} such that 

3 

XjVAjizo) + 71 Vy(0o) + j2VK{zo) = 0. (3.12) 

j=i 

Hence, in general, the location of concentration points might depend also on the 
(fixed) external electromagnetic potential A. If &* = 0, then (Se) does not admit 
any sequence of bound-state solutions concentrating somewhere pointwise. 

Corollary 3.2. The location of energy- concentration points of a sequence of bound- 
state solutions to problem (S^) is independent of the external electromagnetic field 
B (and there holds Xj = for all j = 1, 2, 3 in (3.12) ).If& = 0, then {Se) does not 
admit any sequence of bound-state solutions concentrating somewhere energetically. 

Remark 3.1. Despite the fact that both pointwise and energy concentration are 
gauge invariant, the necessary condition (3.12) is not, in general, unless A^- = for 
all j = 1, 2, 3. Hence, it seems natural to conjecture that the answer to Question 2.1 
is always affirmative. 

Corollary 3.3. Assume that f{u) is such that, for every z G R^, problem (1.3) 

admits a unique positive radial solution, up to translations. Then, if z is an energy- 
concentration point it is a classical critical point ofYjj.. 
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We refer the reader to [5, Theorem 2.5 and Theorem 4.2] for some results ensuring 
uniqueness for (1.3) under some additional hypothesis on f{u). 

We finish the paper with a simple but interesting property of the family {&p}pe{i,5) ■ 
Observation 3.1. Assume that f{u) satisfies (2.6) and that 

limsup-^^^^^^ < DO and liminf |Vi(r(x)| > 0. 

We denote by Crit(-ftr) the set of critical points of K, which is a compact set in 
light of the above assumption. Then, it is a simple task to check that 

lim distR3(6p,Crit(i<')) = 0, 
that is, if p is close to the critical exponent 5, the spikes locate close to Crit(i^). 
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